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ABSTRACT 

We observed 42 molecular condensations within previously identified bright-rimmed clouds 
in the ammonia rotational inversion lines NH3(1,1), (2,2), (3,3) and (4,4) using the Green Bank 
Telescope in Green Bank, West Virginia. Using the relative peaks of the ammonia lines and 
their hyperfine satellites we have determined important parameters of these clouds, including 
rotational temperatures and column densities. 

These observations confirm the presence of dense gas towards IRAS point sources detected 
at submillimetre wavelengths. Derived physical properties allow us to refine the sample of 
bright-rimmed clouds into those likely to be sites of star formation, triggered via the process 
of radiatively-driven implosion. An investigation of the physical properties of our sources 
show that triggered sources are host to greater turbulent velocity dispersions, likely indicative 
of shock motions within the cloud material. These may be attributed to the passage of triggered 
shocks or simply the associati on of outflow activity with the sources. 

In all, we have refined the lSugitani et al.l d 1 99 lb catalogue to 15 clouds which are clearly 
star-forming and influenced by external photoionisation-induced shocks. These sources may 
be said, with high confidence, to represent the best examples of triggering within bright- 
rimmed clouds. 
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1 INTRODUCTION 

Bright-rimmed clouds (BRCs) are small molecular clouds asso- 
ciated with old (>1 Myr) HII regions and are potential exam- 
ples of the radiatively driven implosion (RDI) mode of triggered 
star formation. In the RDI model the UV flux of associated OB 
star(s) powering the HII regions ionizes the external layers of small 
clouds. The ionisation of the cloud's surface leads to the forma- 
tion of an ionized boundary layer (IBL), which subsequently ex- 
pands due to the increased pressure of the ionised gas. The IBL 
expands into the intercloud medium and an ionisation front pre- 
ceded by a shock in the neutral gas prop agates into the cloud 
dLefloch & Lazareffll 1994 1 White et alll997h . The associated pres- 
sure often sweeps the molecular material of the cloud into a 
cometary morphology with a dense core located at the 'head' of 
the cometary globule. 

If the IBL has a pressure greater than, or equal to, the turbu- 
lent and thermal pressure of the cloud then photo-ionisation shocks 
and a D-critical ionisation front propagate into the cloud's interior. 



E-mail: lkm@astro.livjm.ac.uk 



These shocks are driven into the clouds; compressing and heating 
the molecular gas, leading to th e formation of dense cores and po- 
tentia lly causing their collapse l lBertoldilll989l : lLefloch & Lazarefj 
1 1994b . This mechanism could be responsible for triggering the for- 
mation of several hundred stars for each HII region, possibly ac- 
count ing for as much as 10-15% of the low to intermediate mass 
stars dSugitanietalJl99lh . The morphologies and orientations with 
respect to nearby OB stars and optically bright rims identify BRCs 
as potential sites of triggered star formation. Moreover, their rela- 
tive isolation and simple geometry make BRCs an ideal laboratory 
in which to study this mode of triggered star formation. 

Sugitani et al. (1991, 1994) created a catalogue of 89 
BRCs found to be associated with IRAS point sources that have 
IR colours consistent with being protostellar in nature. There 
have been a numbe r of recent studies of sm all sub-samples 
jPe Vries et ai1l2002l; lUrquhart et all [20041. l2006h an d individual 
sources ^Thompson et alj|2004l ; lUrquhart et al J 120070 taken from 
this catalogue. All of these have found evidence consistent with the 
hypothesis that the star formation within the observed sources has 
been triggered; however, this tells us little about the more global 
aspects of this mode of triggered star formation. 
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Figure 1. Multi-wavelength composite image of SFO 13, Spitzer IRAC 
bands at 4.5 and 8.0 /jm are coloured green and blue respectively, a Spitzer 
MIPS 24 /im image is incorporated in red. 21 cm NVSS contours are over- 
laid in white, beginning at a level of two times the r.m.s. noise (cr) of the 
radio image and increasing in units of 1 cr. SCUBA 850 /im contours are 
overlaid in yellow, beginning at a level of 6 cr and increasing in units of 3 cr. 



We have therefore designed a multi-wavelength programme of 
observations to investigate the current level of star formation taking 
place within the whole sample of clouds, and to ascertain whether 
the interaction with the HII region has been a contributing factor 
in that formation. Th e programme includes the use of radio con- 
tinuu m observation s l Thompson et al.l |2004|) and NVSS archival 



data l lMorgan et al.ll2004l) to trace the ionised gas associated with 
the IBL, CO molecular line observations to probe the kinematic s 
of the protostellar cores dUrquhart et al.l2009l ; TMorgan et al.l2009h . 
and submillimetre observations to trace the dust associated with 
the protostellar enve lopes and obtain core temperatures and masses 
jMorgan et alj2008h . In addition to our programme of observations 
we have obtained archival near- to far-infrared imaging and pho- 
tometry (2MASS, MSX, IRAS, Spitzer IRAC and MIPS) in an ef- 
fort to build up a detailed picture of the structure within BRCs and 
a census of the current state of star formation within them. 

A key element of our campaign is to answer the question of 
whether the observed star formation within BRCs has been trig- 
gered by the action of the OB star or is pre-existing and is simply 
being unveiled by the dispersal of the cloud. In order to address 
this question, we must first confirm the existence of protostellar 
cores within BRCs and ch aracterise their proper ties. The sample 
of 44 BRCs catalogued by ISugitani et al.l d 199 lh forms the basis 
for our campaign in the northern hemisphere. A high proportion of 
these clouds have been found to be undergoing a strong interac- 
tion with their surrounding HII region l lMorga n et alj2004|) and ar e 



associated with embedded protostellar cores I Morgan et alj|2008h 



Furthermore, a significant proportion of these sources have been 



found to be associated with water masers fcesaroni et al.l 1 198S 


Henning et al.ll992l;IWouterloot et alj 


19931: IXiang & Turned 1995 


Claussen et al.ll 19961: IValdettaro et al. 


120051, 120081) and kinematic 



line p rofiles indicative of molecular outflows (see iMorgan et al .1 
2009; and references therein). 



In Fig.Q]we present an example of a BRC that appears to be a 
good candidate for triggered star formation. In this figure we have 
combined two Spitzer IRAC bands and a MIPS image to create a 
three colour image of SFO 13 (IRAC 4.5 and 8.0 ,um and MIPS 
24 yum bands are coloured in green, blue and red respectively). We 
have over-plotted contours of the NVSS radio continuum emission 
(white) and the SCUBA 850 j-im emission (yellow) which trace 
ionised gas and thermal dust emission respectively. 

The multiple wavelength data that comprise Fig.[T]show many 
different features of this small molecular cloud on the edge of an 
HII region. T he 4.5 /im emission i n green likely indicates shocked 
H2 gas (e.g. iTeixeira et al .1 12008|) , possibly associated with the 
jets/outflows of young stars. Probing the interface between the 
neutral molecular material of the cloud itself and the ionized gas 
of the HII region, the 8.0 yum emission, traced in blue, indicates 
the presence of polycyclic aromatic hydrocarbons (PAHs), com- 
plex m olecules which absorb far UV photons and re-emit in the IR 
regime dLeger&Pugetlll984l) . Dominant lines of these PAHs are at 
7.7 and 8.6 /im which are within the passband of the Spitzer 8.0 /jm 
band. This reveals the photon-dominated region (PDR), which is 
preceded by the ionization front at the head of the cloud. It is here 
that the cloud is activel y supporting the presence o f an established 
IBL (as described in e.g. lLefloch & Lazareffl ( fl994l) ). emitting free- 
free emission which is traced by the 20 cm NVSS contours in white. 
The 24 /jm emis sion, shown in red in Fig.Q] traces small grains of 
heated dust (e.g. iKoenig et al.l2008h . This emission peaks in a cen- 
trally condensed region at the head of the cloud, coincident with 
the peak of 850 yum contours in yellow, which show submillimetre- 
emission associated with warm dust within the cloud and surround- 
ing the protostellar core at the head of the cloud. 

In order to explore the possible origins of these protostars, it is 
essential to explore the physical properties of the protostellar cores 
themselves and their environs. In this paper we present a set of am- 
monia observations made towards 42 BRCs in an effort to confirm 
the presence of protostellar cores and to probe their physical condi- 
tions. Analysis of the line properties of ammonia rotational transi- 
tions at ~ 23 GHz allows for determinations of kinetic temperature, 
column density and velocity dispersion. These properties will allow 
the definitive identification of protostellar cores within our sample 
and allow us to determine the global properties of molecular mate- 
rial within BRCs. 

The structure of the paper is as follows: in Sect. 2 we describe 
the observational set up, data reduction procedures and derive the 
physical parameters. In Sect. 3 we present the observational results, 
in Sect. 4 we discuss these results with respect to triggered star 
formation. In Sect. 5 we present a summary and highlight our main 
findings. 



2 GBT OBSERVATIONS 
2.1 Observational set up 

Observations of the NH 3 (1,1) and (2,2), (3,3) and (4,4) rotational 
transitions were made over three sessions from the 11''' of February 
to the 20''' of March 2005 on the Green Bank telescope operated by 
the National Radio Astronomy ObservatorjQ. 



1 The National Radio Astronomy Observatory is a facility of the National 
Science Foundation operated under cooperative agreement by Associated 
Universities, Inc. 
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Table 1. Observational setup 

Sate Centre Frequencies (GHz) RMS" (K) T;' ys (K) 

11"' Feb 2005 23.710, 23.877 O05 41.0 

12"' Feb 2005 23.710,23.877 0.04 31.2 

20"' Mar 2005 23.695,23.724,23.871,24.140 0.07 38.4 

"-RMS and Tsys values represent mean values of measurements taken 
from single two minute integrations at intervals over the entire observing 
period. 



Observations were performed in 'nod' mode in order to re- 
move sky contributions and a noise diode was observed in switch- 
ing mode throughout the observations in order to achieve abso- 
lute flux calibration on the scale to an accuracy of ~10%. 
Atmospheric opacity values (at zenith) were determined from 
local weather models, collated in the archive maintained by 
R.Maddalena. Two spectrometer set ups were used during the ob- 
servations due to hardware issues at the time of observation, which 
resulted in the first two sessions being observed with poorer spec- 
tral resolution than the third. Spectrometer problems meant that 
only two spectral windows could be simultaneously observed in 
the first two sessions. The observed bandwidth was thus increased 
so that the (1,1) and (2,2) transitions could be observed in the same 
50 MHz (~650 km s _I ) bandpass with the (3,3) transition in another 
spectral window (the (4,4) transition was not observed in these ses- 
sions). The first two sessions were observed with a spectral channel 
width of 6. 1 kHz (~0.08 km s ). By the time of the third observing 
session, hardware problems had been resolved and all four transi- 
tions were observed at the optimum spectral resolution of ~3 kHz 
(-0.04 km s" 1 ) in 12.5 MHz bandwidth (-160 km s" 1 ) spectral 
windows. All transitions were observed simultaneously using the 
configurations presented in Table [T] 

Weather conditions were stable across all three observing ses- 
sions, typical pointing offsets were -5" in azimuth and eleva- 
tion with a beamsize of -29". Sources were observed for periods 
ranging between four and 66 minutes, dependent upon the signal 
strength of observed lines. The median observing time was 20 min- 
utes per object. 

In total, observations were m ade towards 42 BRCs, in cluding 
SFO HE and 11NE, identified bv lThompson et alj d2004l) but not 
included in the original SFO catalogue. Of the 44 sources listed in 
the northern SFO catalogue, four (SFO 2, 3, 27 and 40) were not 
observed due to time and hour angle limitations. The data of two 
sources (SFO 1 1 and 32) were later found to be corrupt, resulting 
in a total source sample size of 40 out of a potential 46. In Table [2] 
we present a summary of the BRCs observed and the pointing cen- 
tres used for the observations. Source positions were largely deter- 
mined from examination of SCU BA o bserved dust emission m aps, 
presented in lMorgan et al.1 d2008l) and lThompson et all d2004l) . For 
sources not observed in those works, the coordinates of the IRAS 
source associated with the relevant object were used. 

2.2 Data reduction 

Data were reduced using the GBTIDL0 data analysis package. Bad 
scans were removed and channels outside the region of interest 
were discarded. In reduction, use was made of the GBTIDL pro- 
cedural ability to smooth the 'off'. This involves the smoothing of 

2 http://gbtidl.nrao.edu/ 



Table 2. Summary of source names, observed positions and observation 
dates. 



Source 


IRAS 
Source 


a 

(J2000) 




S 

(J2000) 


Obs. Date 


SFO 


01 


23568+6706 


23 


59 


32 


5 


67:24:06 


20/3/2005 


SFO 


04 


00560+6037 


00 


59 


01 


6 


60:53:27 


1 1/2/2005 


SFO 


05 


02252+6120 


02 


29 


01 


6 


61:33:33 


1 1/2/2005 


SFO 


06 


02309+6034 


02 


34 


45 


1 


60:47:49 


20/3/2005 


SFO 


07 


02310+6133 


02 


34 


48 





61:46:36 


1 1/2/2005 


















20/3/2005 


SFO 


08 


02318+6106 


02 


35 


37 


5 


61:19:37 


20/3/2005 


SFO 


09 


02326+6110 


02 


36 


27 


6 


61:24:01 


1 1/2/2005 


SFO 


10 


02443+6012 


02 


48 


12 


2 


60:24:36 


1 1/2/2005 


SFO 


11 


02476+5950 


02 


51 


33 


7 


60:03:54 


20/3/2005 


SFO 


1 IE 




02 


51 


57 


1 


60:06:55 


20/3/2005 


SFO 


11NE 




02 


52 


14 


9 


60:02:54 


20/3/2005 


SFO 


12 


02511+6023 


02 


55 


01 


1 


60:35:45 


12/2/2005 


SFO 


13 


02570+6028 


03 


00 


55 


4 


60:40:18 


1 1/2/2005 


SFO 


14 


02575+6017 


03 


01 


31 


5 


60:29:20 


1 1/2/2005 


















20/3/2005 


SFO 


15 


05202+3309 


05 


23 


28 


3 


33:11:49 


1 1/2/2005 


SFO 


16 


05173-0555 


05 


19 


48 


1 


-05:52:04 


20/3/2005 


SFO 


17 


05286+1203 


05 


31 


28 


1 


12:05:21 


20/3/2005 


SFO 


18 


05417+0907 


05 


44 


29 


g 


09:08:55 


12/2/2005 


















20/3/2005 


SFO 


19 


05320-0300 


05 


34 


30 


6 


-02:58:15 


20/3/2005 


SFO 


20 


05355-0146 


05 


38 


04 


8 


-01:45:10 


20/3/2005 


SFO 


21 


05371-0338 


05 


39 


41 


2 


-03:37:12 


20/3/2005 


SFO 


22 


05359-0515 


05 


38 


26 


2 


-05:14:06 


20/3/2005 


SFO 


23 


06199+2311 


06 


22 


58 


2 


23:10:11 


12/2/2005 


















20/3/2005 


SFO 


24 


06322+0427 


06 


34 


52 


7 


04:25:21 


12/2/2005 


SFO 


25 


06382+1017 


06 


41 


03 


3 


10:15:07 


12/2/2005 


SFO 


26 


07014-1141 


07 


03 


47 


2 


-11:45:46 


20/3/2005 


SFO 


28 


07023-1017 


07 


04 


43 


6 


-10:21:58 


20/3/2005 


SFO 


29 


07025-1204 


07 


04 


52 


5 


-12:09:25 


20/3/2005 


SFO 


30 


18159-1346 


18 


18 


46 


6 


-13:44:26 


12/2/2005 


SFO 


31 


20489+4410 


20 


50 


42 


9 


44:21:56 


12/2/2005 


SFO 


32 


21308+5710 


21 


32 


29 


5 


57:24:33 


20/3/2005 


SFO 


33 


21316+5716 


21 


33 


13 


7 


57:30:03 


12/2/2005 


















20/3/2005 


SFO 


34 


21320+5750 


21 


33 


32 


2 


58:03:34 


12/2/2005 


















20/3/2005 


SFO 


35 


21345+5818 


21 


36 


05 


5 


58:31:38 


20/3/2005 


SFO 


36 


21346+5714 


21 


36 


07 





57:26:40 


12/2/2005 


SFO 


37 


21388+5622 


21 


40 


28 


9 


56:35:54 


12/2/2005 


SFO 


38 


21391+5802 


21 


40 


41 


8 


58:16:12 


12/2/2005 


















20/3/2005 


SFO 


39 


21445+5712 


21 


46 


06 


7 


57:26:38 


12/2/2005 


SFO 


41 


21448+5704 


21 


46 


29 


3 


57:18:41 


20/3/2005 


SFO 


42 


21450+5658 


21 


46 


36 


8 


57:12:25 


20/3/2005 


SFO 


43 


22458+5746 


22 


47 


50 


1 


58:02:47 


12/2/2005 


SFO 


44 


22272+6358A 


22 


28 


51 


5 


64:13:37 


12/2/2005 



the reference spectra before the usual spectral subtraction and cal- 
ibration are performed. This results in higher signal to noise ratios 
for the final calibrated spectrum with no loss of spectral resolution. 
After careful testing of this procedure, the default smoothing ra- 
tio of 16 channels was used for all sources except SFO 31. This 
source showed emission in the reference spectrum which resulted 
in negative emission seen in the final calibrated spectrum. Through 
extensive smoothing of the reference spectrum (by 250 channels) 
this 'off' emission was greatly minimised. Temperature scale cor- 
rections for atmospheric opacity were made using the zenith values 
provided from local weather models. 
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Table 3. Sources not detected in our survey with RMS values of the obser- 
vations. 



Source 


RMS (mK) 


SFO04 


14.8 


SFO06 


9.6 


SFO08 


22.3 


SFO 10 


12.8 


SFO 15 


15.3 


SFO 19 


33.8 


SFO 21 


20.6 


SFO 22 


19.4 


SFO 28 


19.6 



As data from different observing sessions were taken in dif- 
ferent configurations, in certain cases, individual objects may have 
been observed in multiple configurations. Spectra for these ob- 
jects were smoothed to a common resolution before combination. 
All (1,1) and (2,2) spectra were ultimately smoothed to a resolu- 
tion of 7.4 kHz (~0.09 km s -1 ) while (3,3) and (4,4) spectra were 
smoothed to a resolution of 50.2 kHz (due to the expected lower 
signal to noise ratio). A polynomial baseline was subtracted from 
all spectra for normalisation. Several sources were detected in the 
(3,3) and (4,4) transitions. However, analyses of these data con- 
tribute little to the scope of this paper and so the spectra and related 
parameters are presented in appendix [A] 

The corrected antenna temperature T* A , hyperfine linewidth 
Av, and central velocity Vlsr values for each source in the (1,1) 
and (2,2) inversion transitions were determined by fitting with the 
NH3(1,1) METHOD process in the CLASS data analysis pack- 
ag^. In some cases (SFO 09, HE, 17, 26, 35 and 42) the main 
quadrupole of ammonia was detected but the hyperfines were not 
seen at a level required for CLASS to fit them. For these sources, 
and the (3,3) and (4,4) transitions presented in Appendix [A] the 
main quadrupoles were fitted by a single Gaussian using the 'fit- 
gauss' procedure in GBTIDL. 

2.3 Deriving physical parameters 

The optical depth (t) associated with ammonia emission in our 
sources may be determined through the ratio of main to satellite an- 
tenna temperatures of the (1,1) inversion transition dHo & Townej 
Il983h : 

AT* A (J,K) m _ i- e -rU,K)„ 
AT^l K) s ~ 1 - e- aT «*>" 

where m and s subscripts indicate quantities associated with the 
main and satellite quadrupole lines, respectively, and a is the ratio 
of intensities of satellite to main lines (0.278 and 0.221 for inner 
and outer satellites, respectively). Optical depth values were de- 
termined through use of the NH3(1,1) METHOD process in the 
CLASS data analysis package. 

Given the optical depth associated with the (1,1) transition, 
corrected antenna temperatures can be used to calculate the rota- 
tional tom2erature(7^)_associated with the (2,2) and (1,1) transi- 
tions dHo &Towneil 19831) : 



]n{=om Ui _ ^zm (i _ e -^u,)H 
\T„(i,r, [ ^m(i,i) v IS 



where T = E < 2 ~\ « 41.5K. 

For 7k < 7o, a relationship between rotational and kinetic tem- 
perature may be calculated by consideration of (1,1), (2,2), an d 
(2,1) states only l lWalmslev & Ungerechtslll983l ; ISwift et al]2005l) : 



l + ain[l + 0.6exp(-Jf)]' 

although kinetic temperatures calculated by this method associated 
with T T a T may be overestimated. Calculated T k values range 
from 11.9 to 27. 3K, suggesting that this approximation is reason- 
able for these source objects. 

In order to calculate the column density of our sources, it is 
necessary to calculate the excitation temperature of each source. 
This may be found via 

K = r] mi Vt[J(.T e% ) - J(T hg )][l - e T °>\ (4) 

where the quantity J y {T) is defined as > ^ mb [ s me ma i n 

beam efficiency of the GBT (0. 893) and r\i is the filling factor of the 
source in the relevant observation. The filling factor of each source 
is somewhat hard to determine; values of 7" ex derived assuming a 
value of t]{ equal to unity are presented in Table [4] along with val- 
ues of t]{ determined by assuming LTE in our sources and letting 
Te X = 7k. This assumption is supported by th e rough equivalence o f 
T t , r k and the values of 7 dust determined bv lMorgan et alj J2008l) . 
A typical value of a filling f actor in a similar study was determined 
bv lRosolowskv et all d2008[) in their study of dense cores in Perseus 
to be ~0.3. Given the relatively large distances of our sources com- 
pared to Perseus, our values of r] f appear reasonable, ranging from 
0.03 to 0.35 with a mean of 0.12. 

If excitation conditions are homogeneous along the beam and 
all hyperfine lines have the same excitation temperature, then the 
column densities at a g iven (J,K=J) transition can be written as 
dRosolowskv et alj|2008l) 

47r ,/2 v 3 r hv -,-1 

N(l, 1) = 1 - e~n^ Avt, (5) 

Vln2c 3 A L J 

where A is the Einstein spontaneous emission coefficient. Our de- 
termination of the filling factors associated with each source indi- 
cate that the T cx of our sources may be significantly underestimated. 
This is likely due to beam dilution or clumping within the beam of 
our observations, in order to account for this effect we have as- 
sumed LTE in our sources and set 7* cx = T k in our determinations of 
column density. 

Overall column density is obtained from N(l,l) and/or N(2,2) 
and the partition function: 

Z Z(i) Z(i) £f 

where the partition function, 

-h [BJ(J+ 1) + (C-B)J 2 ] 
Zj = (2J + 1) S(J) exp L - ; '—L. ( 7 ) 

kr k 

The values of the rotational constants B and C are 2981 17 and 
186726 MHz respectively and the function S(J) is 2 for J=3,6,9,... 
and 1 for all other J. 

Values for optical depth, T cx and column density are listed for 
each source in Table|4] 



3 http://www.iram.fr/IRAMFR/GILDAS 



4 from 'The Proposer's Guide for the Green Bank Telescope' 
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Table 4. Observationally measured and derived physical source parameters. 



Source 



NH 3 (1,1) 



NH 3 (2,2) 



' A 

(K) 



Av 
(km s" 1 



(K) 



Av Vlsr 
(kms" 1 ) (kms 4 ) 



1 ex 
(K) 



If 



7k Wnh 3 
(K) (10 13 cnT 2 ) 



Tdust 
(K) 



(10 2l cnr 2 ) 



Triggered Sources 


SFO 01 


0.52±0.01 


0.58 


0.21±0.01 


0.78 


-13.5 


1.44±0.184 


2.9±0.1 


0.07 


19.2±0.5 


26.6±1.8 


26 


11.0 


cpn n^ 
ctrKJ Uj 


U.J /+U.U1 


1 9*7 
1 .Z / 


n i Q-i-H n 1 

U. lo+U.Ui 


1 A 8 


-ji.j 


u.zu±u.uz / 


A 9-i.n 9 
4.Z+U.Z 


n 9i 

U.Z1 


Oil 7 


o A-i_n c 
y.o+u.o 


9 1 

Z4 


9^ 7 


ppn 07 
vtl u u / 


S8+0 01 


1 54 


99-1-0 01 
U.ZZ1U.U I 


1 48 


-39.9 


0.64±0.006 


1 S-i-0 


14 


1 Q fij-O 1 


32 0+0 2 


2 1 
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a-Optical depth, temperature and column density are not determinable for these sources as satellite (1,1) qu adrupole hyperfines a re too weak for detection 
and/or the (2,2) main line was not detected. b-Determined from the observations of Morgan et al. 1 2008). 



3 RESULTS 

We detected ammonia (1,1) emission towards 31 of the 40 suc- 
cessfully observed sources. We present spectral plots of both the 
NH 3 (1,1) and (2,2) lines towards all of these in Fig. [5] In these 
plots we show the baseline corrected data (black histogram) and 
the model fits (red line) to the data. The measured properties of the 
main NH? (1,1) and NH? (2,2) lines are given in Table [4] In this 
table we have separated the sample into those clouds considered 
to be good triggered star formation candidates and those in which 
the star formation is unlikely to have been induced. The determi- 
nation of w hich sources are trigg ered candidates and which not is 
described in iMorgan et al.l j2009t ) and is based upon the detection 
of a strong IBL in radio and PDR in IR emission. 

The hyperfine structure of the NH 3 (1,1) transition is clearly 
resolved towards 25 BRCs, as is the main quadrupole of the 
NH3 (2,2) transition, allowing estimates of the optical depth, kine- 
matic and excitation temperatures and column densities to be 
obtained (Several sources also show hyperfine structure in the 
NH 3 (2,2) transition 07, 14, 16, 18, 30, 38, 44). Towards the remain- 
ing six sources only the main NH 3 (1,1) quadrupole line is visible 
above the noise and of these the NH 3 (2,2) transition is only visible 
in three. Our results can thus be separated into three groups, strong 



detections, weak detections and non-detections, each of which con- 
tain 25, 6 and 9 clouds, respectively. Non-detections with respective 
RMS noise levels are listed in Table [3] 

The overall detection rate is 75%, with both transitions being 
strongly detected towards roughly ~60% of the sources. Ammonia 
is a high density tracer needing a crit ical density o f n crit =10 4 cirT 3 
before becoming thermally excited fSwade 1989). Due to the re- 
quirement for this density before ammonia becomes thermally ex- 
cited, ammonia is better able to probe the properties of the dense 
protostellar core surrounding the accreting protostar than CO or 
dust continuum emission, which tend to probe the whole column 
of the gas along the line of sight. 

As previously mentioned, the ammonia observations were 
centred on the positio n of either the dust cores revealed by the 
SCUBA observations jMorgan et al]|2008l) . or towards the IRAS 
point source thought to be protostellar in nature. It is somewhat 
surprising then to find 25% of sources observed resulted in non- 
detections, which would seem to rule out the presence of a dense 
core and subsequently call into question the association of these 
clouds with protostars. Moreover, the weaker emission detected to- 
wards another six clouds may indicate that conditions are not con- 
ducive for star formation. We will investigate these weak and non- 
detections in more detail in Sect. 14. II 
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Figure 2. Histograms showing the distributions of the physical parameters derived from the ammonia detections. Filled red histograms represent our entire 
sample while blue represents our triggered sample only. 



Table S. Summary of measured and derived parameters. 



Parameter 


Min 


Max 


Mean 


Median 


Std. dev. 


Excitation Temp. (K) 


2.5 


4.3 


3.2 


3.1 


0.5 


Kinetic Temp. (K) 


11.9 


27.3 


19.5 


19.9 


4.6 


T 


0.20 


3.55 


1.15 


0.75 


0.98 


Log[N(NH 3 ) (cm" 2 )] 


13.9 


15.2 


14.4 


14.4 


0.3 


Av (km s -1 ) 


0.48 


1.93 


1.07 


1.10 


0.38 



In Fig. [2] we present a set of four histograms showing the dis- 
tribution of ammonia (1,1) linewidth, optical depth, kinetic tem- 
perature and column density. In these plots we show the distribu- 
tion for all of the detections (filled red) and the para meters associ- 
ated w ith BRCs identified as triggered candidates bv lMorgan et al.l 
j2009h (filled blue). The ranges and averages of the properties of 
all detected sources are presented in Table[5] 

The histograms of Fig. [2] show that the temperatures of our 
sources do not show any large variation dependent upon their trig- 



gered status. It would appear that the fundamental observational 
properties of ammonia emission are somewhat insensitive to en- 
vironmental factors. Temperatures and linewidths of protostellar 
sources, as observed in ammonia, appear fairly constant across 
many s ample bases. Some e xamples drawn from the literature 
include iMolinari et ail d 19961) (A^ =1.76 km s _1 , 7\=22.0 K, 260 



IRAS sources), iJiiina et al 



19991) (Av=0.74 km s" 1 , T r =14.7 K, 
264 dense core s drawn from the literature, 1971-1999) and 
IWu et alJ l l2006l) (Av=1.54kms-', 7\=19.0K, 27 IRAS sources 
associated with water masers). 



4 DISCUSSION 

We used the results of our infrared, submillim etre, molecular line 
and radio analyses jMorgan et al.ll2004l2008l , 2009) to identify 26 
BRCs that show strong evidence they are 1) undergoing recently 
initiated star formation and 2) are being subjected to intense lev- 
els of ionizing radiation. These clouds are therefore considered to 
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be excellent candidates in which the observed star formation may 
have been triggered via RDI. Although there is evidence that star 
formation is present in many of the remaining 18 clouds of the SFO 
catalogue there is little or no evidence that the ionisation from the 
nearby OB stars is having a significant dynamical impact on these 
clouds. In this section therefore we will only concern ourselves 
with the BRCs that remain in the refined triggered BRC catalogue 
presented in Morgan et al. (2009). 



4.1 Ammonia detections towards triggered sources 

In Morgan et al. (2009) we identified 26 clouds from the SFO 
catalogue in which there is a strong likelihood that observed star 
formation is the result of the RDI triggering process. Twenty-one 
of these have been observed in ammonia here. Four of the trig- 
gered candidates are non-detections (SFO 4, 6, 10, 15) and a fur- 
ther two fall into the weak category (SFO 35 and 42). These six 
clouds were e ither not detected or only marginally detected by 
iMorgan et all (2008) at submillimetre wavelengths. SFO 6 was 
not detected in that study and the remaining five sources were 
marginally detected at 850 Li.m but not detected at 450 ixm. These 
marginal/non-detections suggest that these clouds are not host to 
protostars, though they may still be at an early point in their evolu- 
tion with respect to the observed ionization fronts and may, in the 
future, develop protostellar cores. These sources shall be discarded 
from further consideration of triggered protostellar souces. 

Ten of the sources in our observational sample have been pre- 
viously associated with water maser activity (see Table [6j. These 
account for 40% of our 'strongly detected' sources, a figure con- 
sis tent with the 40% detection rate for Class protostars found 
bv lFuruva et alj J200lb . All ten of the sources associated with wa- 
ter masers are also associated with outflow activity (see Table [6), 
supporting the general associat ion of jets and outflows with wa- 
ter masers dCodella etai]|2004l and references therein). Of these 
ten sources, only two are found within our non-triggered sam- 
ple. While this finding is perhaps not statistically significant, given 
the small sample size, it does indicate that high radiation environ- 
ments are unlik e ly to suppress maser emission, as suggested by 
Valdettaro et al. | d2008h . The presence of molecular outflows and 
water masers towards so many triggered candidates is a strong in- 
dication that star formation is currently taking place within them. 
Observed linewidths within sources host to outflows and masers 
are typically greater than in other sources (median linewidths are 
1.4 km s" 1 for maser sources, compared to 1.1 km s -1 for other 
sources), indicating that the high density gas surrounding the pro- 
tostellar cores traced by ammonia emission is dynamically linked 
to the collisional processes associated with the maser emission. 



4.2 Comparing protostellar properties derived from different 
tracers 

In this section we will compare the derived ammonia and dust prop- 
erties to identify correlations and anti-correlations in the data and 
to check for inconsistencies between the different tracers. In Sec- 
tion |2.3| we derived the kinetic temperatures and NH 3 column den- 
sities which are the most readily available quantities we have at 
hand to compare with quantities derived from the dust emission. In 
the next section we will describe how the dust temperatures were 
determined and estimate the Hi column densities. 



Table 6. Outflow and water maser activity in our sources 



Source Av (km s 1 ) Outflow Maser 



SFO 05 


1.27 


1 


3 


SFO 07 


1.54 


1 


4 


SFO 14 


1.62 


2 


5 


SFO 18* 


0.64 


1 


6 


SFO 30 


1.43 


1 


7 


SFO 31 


1.13 


1 


8 


SFO 36 


1.15 


1 


9 


SFO 37 


0.83 


2 


9 


SFO 38 


1.86 


2 


7 


SFO 39* 


1.11 


2 


9 



* BRCs not considered to be triggered candidates. 

References (NB, o nl y most recent reference cit ed): lj Morgan et al.U2009h . 
2jWu et all hoo4). 3jXiang & Turner! j 19951). 4fWouterloot et alj j 1993b. 5 - 
iHenning et alj I 992h. 6 faaussen et all ll996t). jJValdettaro et alj feOOSh . 
8 jCesaroni et aljjl988l) . 9 jValdettaro et all J2008b 

4.2.7 Dust properties 

In an earlier paper dMorgan et alll2008h we presented observations 
of submillimetre emission which tracing the distribution of warm 
dust. Spectral energy distributions (SEDs) were determined by fit- 
ting greybody functions to the measured submillimetre fluxes and 
mid- and far-infrared fluxes. Fits were presented for all sources for 
which good quality data was available. 

The H 2 column density associated with each SFO object de- 
tected bv lMorgan et alj j2008l) may be calculated using 

N(H 2 ) = S v /[£lfim H K v B v (T d )] (8) 

where S v is the 850 [im flux density, Si is the solid an gle associated 
with t he 30" aperture used to observe each core in IMorgan et alj 
d2008l) . fx = 2.3 is the mean molecular weight, m H is the mass of a 
hydrogen atom, k v i s the dust opacity per unit mass at 850 iim (0.02 
cm- g~ ', following IMorgan et alj J2008h ) and B v (r d ) is the Planck 
function, evaluated at dust temperature Resulting values of IT 
column density are presented in the final column of Table [4] 

4.2.2 Dust and gas temperatures and column densities 

In Fig. [3] we present two scatter plots comparing the temperatures 
(upper panel) and column densities (lower panel) derived from the 
two tracers of ammonia and submillimetre emission. We show the 
linear-square fit to the data as a solid line. Additionally, in the tem- 
perature plot, we include a dashed line indicating the position of 
the data if both temperatures were equal. There is quite a lot of 
scatter in the distributions seen in both plots. However, there is a 
general correlation between the kinetic and dust temperatures and 
the H 2 and NH 3 column densities. The dust temperatures are gen- 
erally slightly higher than the observed kinetic temperatures. The 
ratio of dust temperature to 7\ ranges from 0.9 to 1.6 with a mean 
of 1.2, this slightly higher dust temperature may be attributed to the 
fact that submillimetre emission is associated with a wide range of 
densities, covering the star-forming core itself as well as the warm 
envelope surrounding the core (and the interface between the two). 
As ammonia emission requires a critical density of ~10 4 cm" 3 , it 
is likely to trace the inner, more dense regions of the protostellar 
core. It should be noted that the dust temperatures were derived 
using fluxes from the IRAS with a significantly larger beam than 
the present observations. These observations are therefore likely to 



8 L. K. Morgan et al. 



\? 


30 


CD 




71 


25 


O 




o: 




Q_ 


20 




0) 

h- 




o 


15 


<D 




C 






10 




the variation of fractional abundance from source to source. 
Overall, the determined values of fractional abundance are typical 
across a wide range of protostellar environments, fro m low-mass 
starless cores jTafalla et alj|2006l: ICrapsi et alj|2007h and low to 
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Figure 3. Upper panel: Dust temperature vs. kinetic temperature, a line of 
best fit is shown as a solid line while the dotted line indicates the line of 
equality. Lower panel: NH3 vs. H2 column density, a line of best fit is shown 
as a solid line while the dotted line indicates the line of constant fractional 
abundance. 



incorporate more of the hotter dust at the edges of the BRCs. We 
therefore conclude that the two tracers are probing material at sim- 
ilar temperatures (difference between median averages is <3 K), 
though the submillimetre observations may incorporate some addi- 
tional material associated with the warm protostellar envelopes and 
cloud rims. 

Although there is some correlation seen in the comparison 
plot of the H 2 and NH 3 column densities the scatter is significant. 
The fractional abundance of NH3 may be found through simple 
comparison of the derived H2 and NH3 column densities. Looking 
at individual sources we find the fractional abundances range 
from a few times 10" 9 to a few times 1(T 8 . The mean fractional 
abundance of NH3 to H2 is 2.6 x 1CT 8 , this is the value used 
in Fig. [3] to illustrate a line of constant fractional abundance. 
The scatter in the plot of NH3 vs. H2 column density reflects 



Note that the beam size of the SCUBA observations was ~ 15", sig- 
nificantly smaller than the ~ 29" beamsize of the ammonia observ ations 
presented here. However, the fluxes taken from lMorgan et all j2008l) were 
summed over a 30" aperture, thus making a direct comparison valid. The 
disparity between the pointing positions of each set of observations is typi- 
cally small (~ 1-2" ) and only exceeds a half-beamwidth for one source, SFO 
29. 



interm ediate mass dense cores dffotzel et alT 2001 : Friesen et al.l 
120091) . to complex, PDR-associate d regions |Larsson et alj|2003l) 



and h igh-mass star forming regions dKuiper et aljl99 Pillai et ahl 
120061) . 

The fractional abundances found here reflect a more general 
trend in the properties of ammonia in star forming regions. The 
physical properties of our sources, as determined from our ammo- 
nia observations, are typical in most star forming environments, 
with only very hot cores showi ng any significant variation in col- 
umn density or temperature (c.f. lLongmore et al1l2007l ; |Pillai et al.l 
120071) . The implication of our analysis is that ammonia, once ex- 
cited beyond its critical threshold, is insensitive to environmen- 
tal circumstances, i.e. resistant to depletion in cold, dense cores 
and likely shielded from photoionisation in high-radiation environ- 
ments. 



4.2.3 Turbulence in BRCs 

The contribution of turbulent motions to the observed linewidth 
of each source was estimated by removing linewidth contributions 
from other broadening agents. Thermal contributions to each line 
were calculated via dVtherm = -^8k B 7\ln2/m NHl , where T k is the 
kinetic temperature of the source in question and iriNH, is the mean 
molecular mass of an ammonia molecule (17.03 amu). 

After removing the thermal broadening contribution to our ob- 
served linewidths we may evaluate the contribution of turbulent 
motions in our clouds. An analysis of the turbulent velocity dis- 
persions (cr = Av > 2 /(81n2)) of our sources, separated based 
upon their triggered status (see Sect.[T]l, reveals some interesting 
differences between the two samples. A histogram of the velocity 
dispersions of the potentially triggered and non-triggered samples 
shows that our non-triggered sources have typically lower velocity 
dispersions than our triggered candidates (Fig. [4}. A Kolmogorov- 
Smirnov test of the cr of each of the two samples indicates that 
the two distributions are drawn from separate populations with a 
probability of 99.8%. An interesting point to note in the distribu- 
tions of the tr of each sample is that the non-triggered sources are 
largely subsonic, based upo n an estimate of the so und speed in our 
sources of ~0.3 km s~' (e.g. lThompson et al.l2004) . In contrast, the 
triggered sample typically exhibit supersonic velocities. 

Several explanations of the difference in cr between our two 
samples exist, the sources in our triggered sample are likely un- 
dergoing the progression of shocks through the host clouds. As 
this sample was selected on the basis of the presence of an ion- 
isation front, the presence of shocks in the observed material 
w ould thus be con s istent with the mode l s of tr i ggering put f orwar d 
bvlBertoldil Jl989h : lLefloch & Lazarefj Jl994h ; lMiao et all J2006h : 
iGritschneder et alj j2009h : lMiao et al.N2009n . An alternative expla- 
nation would be that the increased turbulent motions within our 
triggered sample are due to increased systematic motions within the 
relevant clouds (due to larger numbers of outflows, for example). 
Alternatively, these clouds may simply be larger in extent and/or 
mass. However, we found no trend toward higher dust mass for 
these clouds from the submillimetre observations of iMorgan et al.l 
d2008h . The source luminosities from those observations show that 
the most luminous sources of the enti re sample are those of the 
triggered sub-sample, as suggested by IMorgan et af] J2008h . The 



Ammonia observations of bright-rimmed clouds 9 



Figure 4. Histogram of turbulent velocity dispersion for triggered and non- 
triggered Sources, Group I (triggered) sources are shown in blue and Group 
II sources are shown in red. Binsize is 0.1 km s -1 . 
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median luminosity of the triggered sample is 201 L G while the non- 
triggered sample median average is almost an order of magnitude 
lower at 28 L Q . While this may go some way to explain the higher 
velocity dispersions within the triggered sample, no direct correla- 
tion between source luminosity and linewidth can be drawn. Such 
a correlation might be expected if the higher turbulent velocity dis- 
persi ons of the triggered sample were due to outfl ow momentum 
(e.g. ICabrit & Bertout|[l992l ; iBontemps et alJI 19961) . The ambigu- 
ity in the true nature of the increased velocity dispersions in these 
clouds cannot be resolved with single pointing data and must be ad- 
dressed with maps of the molecular emission associated with each 
region. 



5 SUMMARY AND CONCLUSIONS 

We present observations made with the GBT of the ammonia (1,1) 
and (2,2) inversion transitions tow ards 40 bright-rimme d clouds 
taken from the sample compiled bv lSugitani et all dl99ll) . We de- 
tected the (1,1) transition towards 31 of the sources observed, with 
the (2,2) transitions seen towards 28 of these. We use these emis- 
sion lines to derive the optical depths, kinetic temperatures and NHj 
column densities towards the embedded protostars located within 
these clouds. 

Across the entire sample, the ammonia kinetic temperatures 
of our sources correlate reason ably well with t he dus t temperatures 
seen towards the same cores bv lMorgan et all J2008b . Dust temper- 
atures are typically slightly higher than the observed kinetic tem- 
peratures, likely due to the slightly different material traced by the 
two sets of observations. A correlation between H2 and NH3 col- 
umn densities is also seen, although there is somewhat more scatter 
in the plotted points. This likely reflects variation in the fractional 
abundance of ammonia from source to source. A mean fractional 
abundance of NH3 to H2 is 2.6 x 10 s , typical for most observa- 
tions of protostellar regions. 

Using a combination of mid-infrared, submillimetre and ra- 
dio images, and CO molecular lin e data presented in earlier works 
dMorgan et ai] |2004 120081, |2009|) . we have refined the original 
sample of 44 BRCs. Our efforts have identified 26 bright-rimmed 
clouds in which the data are consistent with the hypothesis that any 
observed star formation is likely to have been triggered. Within this 
refined sample we have detected strong ammonia emission towards 
15. In combination with the submillimetre continuum and CO line 



emission results o f iMorgan et all d2008l) and lMorgan et alj J2009I) 
respectively, in addition to outflow and maser detections from the 
literature, our ammonia detections leave little doubt of the star- 
forming nature of these sources. Having already established the 
likelihood that these sources represent photoionisation-triggering 
processes in progress, these 15 sources are some of the best exam- 
ples yet known of the RDI process. 

An investigation of our samples, separated based upon trig- 
gered status, indicates a bimodality within observations of the tur- 
bulent velocity dispersion. Those sources which have been identi- 
fied as likely triggered in nature show typically supersonic turbulent 
velocity dispersions. While non-triggered sources are more often 
found to be subsonic. This disparity may stem from the presence of 
shocks, traversing the clouds in the triggered sample. It is also pos- 
sible that the higher observed velocity dispersion in the triggered 
sample is simply due to the higher occurence of outflows found in 
that sample. It is tempting to draw conclusions of differing physical 
processes occuring within our triggered and non-triggered samples. 
However, we are not able to determine the true cause of this finding 
without mapping of molecular emission in these clouds. 
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Table Al. Observationally measured parameters of (3,3) and (4,4) detec- 
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Av 


Source 


(K) 


(kms- 1 ) 


(3,3) Sources 


SFO01 


0.10±0.006 


1.45 


SFO 07 


0.10±0.004 


2.04 


SFO 11NE 


0.07±0.005 


1.55 


SFO 14 


0.38±0.007 


2.24 


SFO 18 


0.03±0.004 


2.20 


SFO 33 


0.05±0.003 


1.45 


SFO 38 


0.17±0.005 


3.03 


SFO 41 


0.04±0.005 


1.73 


(4,4) Sources 


SFO 14 


0.04±0.008 


2.56 
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APPENDIX A: (3,3) AND (4,4) DETECTIONS 

Several sources were detected in the (3,3) and (4,4) rotational 
transitions of ammonia. Spectra of these sources are presented in 
Figs. IA1 1 and IA2I and relevant parameters are listed in Table [ATI 
RMS values of these observations are presented in Table lA2l 

Those sources detected in the (3,3) line of ammonia were ex- 
amined for signs of masing. If masing is occuring in this transi- 
tion then the brightness temperature of that source would be larger 
than the brightness temperature in the (1,1) transition and com - 
parable to the (1,1) kinetic temperature (e.g. iKuiperetalll 19951) . 
However, none of our sources show such symptoms. The (4,4) 
line of ammonia may often be observed toward warmer cores (e.g. 
lLongmore et al.ll2007l ). only one of our sources provided a detec- 
tion in this line, SFO 14. This core is indeed one of the warmer 
sources in our sample at 7\ ~26K. 

This paper has been typeset from a TpX/ ETpX file prepared by the 
author. 
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Figure 5. Spectral lines are presented with corrected antenna temperature, 7"^, plotted against doppler shifted velocity, Vlsr- Multiple lines are plotted on the 
same axis range with the (1,1) transition spectra on the bottom and (2,2) spectra on top. Fitted profiles are overlaid in red, determien. 
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Figure 5. (cont.) Spectral lines are presented with corrected antenna temperature, T^, plotted against doppler shifted velocity, Vlsr- Multiple lines are plotted 
on the same axis range with the (1,1) transition spectra on the bottom and (2,2) spectra on top. 
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Figure Al. NH3 (3,3) transition spectral lines. 7^, plotted against doppler shifted velocity, Vlsr, with GBTIDL fitted Gaussian profiles overlaid in red. 
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SFO 14 




Figure A2. NH3 (4,4) transition spectral line observed towards SFO 14. T^, plotted against doppler shifted velocity, Vlsr> with GBTIDL fitted Gaussian 
profiles overlaid in red. 



Table A2. RMS values of our (3,3) and (4,4) transition observations. 



Source (3,3) RMS (mK) (4,4) RMS (mK) 



SFO 01 


13.3 


4.7 


SFO 04 


13.1 




SFO 05 


12.0 




SFO 06 


4.4 


4.2 


SFO 07 


14.8 


4.6 


SFO 08 


10.1 


9,3 


SFO 09 


14.3 




SFO 10 


12.6 




SFO HE 


10.1 


9.9 


SFO 11NE 


10.3 


5.3 


SFO 12 


4.4 




SFO 13 


11.4 




SFO 15 


13.7 




SFO 16 


4.4 


4.7 


SFO 17 


11.0 


11.7 


SFO 18 


5.5 


4.0 


SFO 19 


16.0 


14.1 


SFO 20 


7.3 


6.6 


SFO 21 


9.6 


8.5 


SFO 22 


9.7 


8.7 


SFO 23 


4.5 


3.9 


SFO 24 


7.2 




SFO 25 


5.7 




SFO 26 


10.0 


9.4 


SFO 28 


7.9 


8.2 


SFO 29 


7.3 


6.2 


SFO 30 


7.9 




SFO 31 


3.8 




SFO 33 


6.7 


5.0 


SFO 34 


6.4 


5.2 


SFO 35 


9.4 


9.8 


SFO 36 


4.5 




SFO 37 


5.2 




SFO 38 


28.0 


5.2 


SFO 39 


4.4 




SFO 41 


7.7 


4.6 


SFO 42 


13.4 


10.5 


SFO 43 


4.8 




SFO 44 


5.4 





■ ■ ■ Source not observed in this transition. 



